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I. INTRODUCTION 



The study of low-energy QCD remains a long-standing and exciting challenge since its 
conception in the early 1970's. Over the years, various models and tools have been devel- 
oped in an attempt to study the non-perturbative regime of the theory of strong interactions. 
Despite the success some of these ideas have had, an intimate understanding of low-energy 
QCD and the hadron spectra remains elusive. Recently, progress toward this end is being 
made at many fronts, one of which is lattice QCD and partially quenched chiral perturbation 
theory (PQ%PT). Lattice QCD is a first principles method to calculate QCD observables us- 
ing numerical techniques. While ongoing progress in this area is impressive, one is presently 
restricted to using light quark masses on the order of the strange quark mass, significantly 
heavier than those in nature. Hence one needs to understand the quark mass dependence 
of QCD observables in order to compare lattice simulations with physical QCD. For suffi- 
ciently small quark masses, chiral perturbation theory (%PT) can be used to extrapolate 
unquenched lattice simulations from lattice quark masses to physical quark masses. How- 
ever, for lattice masses that are presently being used the convergence of the chiral expansion 
is uncertain. Baryon masses as well as other observables have been computed using QCD 
and quenched QCD (QQCD) 0,333- 

In QQCD the determinant arising from the path 
integral over the quark fields is set to a constant, which greatly reduces the computing time 
of lattice simulations. There has been some interest in using chiral effective theories to 
extrapolate these lattice data to lighter quark masses, for example (H El ' but there is no 
rigorous limit of QQCD in which one recovers QCD. 

Partially quenched QCD (PQQCD) has been formulated as an alternative to QCD and 
QQCD US [HI U|. In PQQCD, the valence quarks, which are coupled to external states, 



and the sea quarks, which contribute to the quark determinant, are treated as independent 
fields. Thus the valence and sea quark masses can be different. This freedom is generally 
exploited to make the valence quark masses much smaller than the sea quark masses. The 
low energy effective theory of PQQCD is PQxPT 0, The coefficients of operators 
appearing in the theory are the low-energy constants (LEC's) which characterize the short 
distance physics, while the long distance behavior is characterized by meson loops. The 
LEC's of xPT are all contained in PQ%PT which allows rigorous predictions for QCD using 
PQxPT. Moreover, the ability to independently vary the valence and sea quark masses 
greatly increases the allowed parameter space available to PQQCD. This enables a better 
extraction of the LEC's appearing in PQxPT, as they can be fit to a much larger set of PQ 
lattice calculations. 

These last few years have seen much work and progress in calculating properties of baryons 
in PQxPT 13 0, |H, H H H H H Q. In this work we calculate the masses of the 
octet baryons to next-to-next-to leading order (NNLO) in %PT. We also provide the first 
NNLO calculation of the octet baryon masses in PQxPT. These calculations are performed 
in the isospin limit of SU(3) and we keep the decuplet baryons as dynamical fields. In the 
partially quenched calculation we use three valence, three ghost and three sea quarks, with 
two of the sea quarks degenerate. 
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II. HEAVY BARYON CHIRAL PERTURBATION THEORY 



A. Pseudo-Goldstone Bosons 



For massless quarks, QCD exhibits a chiral symmetry, SU (3)l<S>SU (3) r®U{1)v, which is 
spontaneously broken down to SU(3)v <8> U (l)y. Chiral perturbation theory, the low-energy 
effective theory, emerges by expanding about the physical vacuum state of QCD. In the limit 
of massless quarks, the pions, kaons and eta emerge as the Goldstone bosons of the broken 
SU(3)a- Given that the quark masses are small compared to the scale of chiral symmetry 
breaking, the lowest lying mesons emerge as an octet of pseudo-Goldstone bosons. 

The pseudo-Goldstone bosons are collected in an exponential matrix 



exp 



2 i $ 



7T 



V 



K 



With the above convention, the pion decay constant, /, is 132 MeV. Under an 
SU(3)l <S> SU(3)r transformation these fields transform in the following way 



£ -> L£i? f , £ -> L£tf = U£R\ (2) 

where U is implicitly defined by equations JI} and (j2J). The effective Lagrangian describing 
the strong dynamics of these mesons at leading order in x?T is 

f 2 

£ = ^-Tr (^E^E) + ATr (mj S f + m q Y) . (3) 

Leading order (LO) in the power counting is 0{m q ) and thus 0{m q ) ~ 0(q 2 ), where q is a 
pion momentum. The quark mass matrix, m q , is given in the isospin limit, (m u = = fn), 
by 

m q = diag(m, m, m s ). (4) 
To the meson masses to LO are given by 

2 8A 2 4A . . 2 , , , . 

m^ = —m , m K = — (m + m s ) , = — (m + 2m s ) . (5) 



B. Baryons 

To systematically include the octet and decuplet baryons into the chiral Lagrangian, we 
use heavy baryon (HB^PT) |23i I24I l25i l2fi| , where the baryon fields are redefined in 
terms of velocity dependent fields. 

B v (x) = ^-^-e mBV - x B(x), (6) 

where is the four-velocity of the baryon, B. This field redefinition corresponds to param- 
eterizing the momentum of a nearly on-shell baryon as 

= M B Vf, + kp, (7) 
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where k^ is the residual momentum. The effect of this is to eliminate the standard Dirac 
mass term for the baryons: 



B(i$-M B )B = B v i$B v + 0(—). (8) 

M B 

From Eq. it is easy to verify that derivatives acting on B v bring down powers of k, rather 
than p. Thus, higher dimension operators of the heavy baryon field B v are suppressed 
by powers of Mb- Heavy baryon ;\;PT is applicable in the limit that pion momenta and 
the off-shellness of the baryon are small compared to the chiral symmetry breaking scale, 
A x ~ lGeV, i.e.: 

q, k ■ v < A x . (9) 

In this limit, HB^PT has a consistent derivative expansion as derivatives are suppressed by 
powers of A x . 

When the spin-| decuplet baryons, T, are included in the theory, an additional mass 
parameter, A = Mt — Mb, must be included. The meson masses, m^, where is a pion, 
kaon or eta, are much smaller than A close to the chiral limit, and therefore the decuplet 
resonances can be integrated out of the theory, leaving only the pseudo-Goldstone mesons 
and the octet baryons. Decuplet effects would then show up in the theory as higher dimension 
operators suppressed by powers of C 2 / A (where C ~ 1.5 is the TBn coupling). In the real 
world, A ~ 300 MeV and m n /A ~ 1/2, so 0{m 7T /A) corrections could be sizeable and 
spoil the power counting Additionally, if one is to use present day lattice calculations to 
fit the free parameters of the theory, m l ® tt / A > 1. The situation becomes worse when 
considering (9(m^/A) and 0(m v /A) corrections. There is also ample phenomenological 
evidence in the nucleon sector that suggests the importance of including the lowest lying 
spin- 1 resonances 

To include the decuplet baryons, one must include A in the power counting. The mass 
splitting, A ~ so A and are O(q), where q is a typical small pion momentum, and 
the quark mass, m q , is treated as 0(q 2 ). The baryon mass is treated as M B ~ A x . As 
the LEC's are a priori unknown, we can combine the 1 /Mb and 1/A X expansions into one 
expansion in powers of 1/A X . There is one exception: constraints from reparameterization 
invariance determine the coefficients of some of the higher dimension operators arising in 
the 1/Mb expansion. |3oL l3~D]. Thus these 1/Mb corrections must be kept distinct to insure 
the Lorentz invariance of the heavy baryon theory to a given order. 

The effective Lagrangian at leading order (LO) in the 1/Mb expansion is well known 
In this paper we will embed the baryons in rank-3 flavor tensors. The convenience of this 
choice will become readily apparent when we generalize to PQxPT. The SU(3) matrix of 
the lowest-lying spin-| baryon fields is 

B=| E" 7^-75^° " I ( 10 ) 



2° -7E A ' 



We then embed them in the tensor 



B ijk = ( £ijl&k + tiklRj ) , (11) 
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which has the symmetry properties 

B^k = Bikj and B^k + B^ + B k ji = 0. (12) 

Under chiral 577(3) transformations, this baryon tensor transforms as 

B^k — ► Ui l Uj m U k n Bi mn . (13) 

The spin-| decuplet baryons can be described by a Rarita-Schwinger field, (T^)^, which 
is totally symmetric under the interchange of flavor indices, and which satisfies the con- 
straint, 7 M = 0. We employ the normalization convention that T 111 = A ++ . Under a 
chiral SU(3) transformation, the decuplet tensor transforms identically to the octet tensor: 

T^^U^U^U^T^. (14) 

The octet and decuplet baryon tensors obey the spin-algebra laid out in [23j, which can 
be used to eliminate the Dirac structure of the theory; all Lorentz tensors made from spinors 
can be written in terms of v M and S^, where v M is the four- velocity of the baryon and is 
the covariant spin- vector. 

Thus the Lagrangian to leading order in the 1/Mb expansion can be written as 1 

Z LO = (B iv-D B) + 2a M (BBM+) + 2(3 M (BM+B) + 2a M (~BB) Tr(M+) 
-(T[ivD-A]T^) + 2 7w pA4 + rJ - 2a M (T%)Tr(M + ) 
+2a(BS t *BA ll ) + 2/3 (BS"A„B) + 2H (TS^A^) 

+Vf C W V ^B) + (BA V T»)] . (15) 

Above, Dp is the chiral- covariant derivative which acts on the B and T fields as 

(D»B) ijk = d»B ijk + (V^ l B ljk + (V% l B ilk + (V^ k l B ijt . (16) 
The vector and axial-vector meson fields appearing in the Lagrangian are given by 

V, = \{^ + , A,= l -{id^ - t*d,£), (17) 

and 

m + = - (e t m 9 e t + e^e) • (is) 

In Eq. (|TKj) the brackets, ( ), denote a contraction of the flavor and Lorentz indices. For a 
matrix acting in spin-space and a matrix acting in flavor space, the contractions are 
defined as 

(btb) = B a ' kJl B pm , (r r T„ ) 

(BTY B)= B*' kji Yi Y % 1 B pm , (T T Y Tj 
(BT BY)= B a ' k3i Ti Y k l B pm , (~B T Y» Tj 



For brevity, we will drop the subscript v from the velocity dependent heavy baryon fields. 



-I J- a J- tf,ijk 

- 1 1 a Y i 1 UPW 

= B a ^ l Ti (Y% 1 T mjk . (19) 



5 



Such contractions ensure the proper transformations of the field bilinears under chiral trans- 
formations. To compare with the coefficients used in the standard two-index baryon formu- 
lation |23|, 24[ 26[, it is straightforward to show that: 



a 



Id 

3 



2F 



—D + F, 
3 ' 



% = ~b D + 2b F , f3 M = + b F 



b F + a . 



(20) 



C. Higher Dimensional Operators 

The Lagrangian in Eq. (j!5|) contains some, but not all, terms of 0(q 2 ). To calculate the 
baryon masses to 0(m 2 ) we must include all relevant 0(q 2 ) operators, which contribute to 
the mass via loops, and all relevant 0(q 4 ) operators which contribute at tree level. The 
Lagrangian also includes operators of 0(q 3 ), but they do not contribute to the baryon self- 
energy. 

The higher dimensional operators introduce new LEC's that, in principle, must be de- 
termined from experiments or lattice QCD calculations. However, many of them can be 
derived exactly using reparameterisation invariance (RI) [3(j. The baryon momentum pa- 
rameterization in terms of t> M and /c M is not unique when considering the 1 /Mb corrections. 
When the velocity and residual momentum are simultaneously transformed in the following 
way 

u ->« + -!_ , k^k-e, (21) 
Mb 

the momentum in Eq. ((7j) is unchanged. Reparameterization invariance requires the 
effective Lagrangian to be invariant under such a transformation, which ensures the theory 
is Lorentz invariant to a given order in 1 /Mb- Furthermore, utilizing RI has non-trivial 
consequences as it connects operators at different orders in the 1 /Mb expansion, exactly 
fixing the coefficients of some of the higher dimensional operators with respect to the lower 
ones. We find that the fixed coefficient Lagrangian is 2 



{ 2M B \ Mb J \ M B 



\ Bv ■ A ^ B + ' ( B ^r ■ AB + r^r *J (22) 



where 3 D\ = D 2 - (v ■ D) 2 . 



2 While there are at least three other methods of generating the 1/Mg Lagrangian, the technique of RI 
allows an explicit check for redundant or missing operators. 

3 Using D 2 , instead of D 2 is a choice as one could use the LO equations of motion to eliminate the (v ■ D) 2 
term |33|, (suitably shifting the LECs as well so that observables remain unchanged). This is an aesthetic 
choice as it gives the familiar form of the first relativistic correction to the propagator, k 2 /2Mg. 
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The Lagrangian contains additional 0(q 2 ) operators as well as 0(q A ) operators which are 
invariant under the SU(3) chiral transformations. The operators relevant to the self energy 
of the octet baryons are 

£S 4) = | b? B kjl (A ■ A\ n B njk + b A B kjl (A ■ A) k n B ijn 

+ b A B kj \A^) i l (A% n B lnk + b A B kji B ijk Tr(A-A) 

+ bf B k]l (v ■ A v ■ A), n B njh + bf B kjl (v-AvA) k n B ijn 

+ bl A B hji {v ■ A^v ■ A) 3 n B lnk + bf B hri B tjk Tr(vAvA) 

+ bf B kjt {M+M+)? B njk + 6f B kji (M + M + ) k n B ijn 

+ bfB kj \M + ) l \M + ) J n B lnk + bf B kji B ijk Tr(M + M+) 

+ bfB kji (M+) i n B njk Tr(M + ) + bfB kj \M + ) k n B ijn Tr(M + ) 

+ 6f B kji B ijk Ti (M+) Tr (M+) J . (23) 
The operators of the form (Bv ■ Av ■ AB^j have not been kep t explicit in the literature 



so far for calculations of the octet baryon masses to 0(m ) [2J, |26|, |34j. These operators 



have identical flavor structure to the corresponding (BA ■ AB^ operators. However, these 
two sets of operators have different Lorentz structure, which gives rise to different finite 
mi contributions to the octet baryon self energy calculations. One can choose a renormal- 
ization scheme such that the contributions to the octet baryon masses from these different 
operators can not be distinguished. Therefore, with a suitable redefinition of the b A _ 4 and 
bf_ 7 coefficients, the operators with coefficients b\ A A can be neglected in the baryon mass 
calculations, as their contributions to the masses can be absorbed by the other operators to 
this order in the chiral expansion. 

However, the operators with coefficients b A _ 4 and b\ A A can be distinguished in nN — > ttN 
scattering at tree level, for example. It is therefore useful to keep both types of operators 
explicitly in the Lagrangian, allowing them to be distinguished in the octet baryon mass 
calculation. This also provides a consistent means to determine these LECs, as they are the 
same coefficients which appear in the Lagrangian used for the calculation of other observables 
like ttN — > ttN scattering. 4 

In principle, additional 1 /Mg operators with the same chiral symmetry properties as those 
contained in Eq. (|23|) can be generated. However, these 1/Mg operators do not have their 
coefficients constrained by RI, therefore they can be absorbed by a re-definition of the b A,vA ' M 
coefficients. For example, the Rarita-Schwinger field used to describe the decuplet baryons 
contains un-physical spin-1/2 degrees of freedom which can propagate when the decuplet 
baryons are off mass-shell. It can be shown that the contribution to the octet baryon masses 
from these un-physical degrees of freedom is suppressed by 1 /Mb Therefore these 

effects are implicitly included in the operators with coefficients bf^ A . 



One could also add operators of the form (BS ■ AS ■ AB) , but it is straight forward to show that they are 
a linear combination of the operators in Eq. I|23fl . and therefore not distinct. 
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FIG. 1: One loop graphs which give contributions to M B . The single, double and dashed lines 
correspond to octet baryons, decuplet baryons and mesons, respectively. The filled squares denote 
the axial coupling given in Eq. I|15|). 

The full Lagrangian relevant to the calculation of the octet baryon masses to 0{m 2 q ) is 

Zmi = 2lo + + ££ 4) • (24) 



III. BARYON MASSES 



The mass of the i th baryon in the chiral expansion is 



M, 



B, 



M (p) - MjJ (//) 



M 



(3/2) 
Bi 



(/i) - Mg J (//) + ... 



(25) 



Here, M (//) is the renormalized mass of the octet baryons in the chiral limit, is independent 
of m q and also B^. Mg} is the contribution to the i th octet baryon of the order mf, and /i 
is the dimensional regularization scale. For this calculation we use a renormalization scheme 
of always subtracting terms proportional to 



1 



— 7 + 1 + log4-7T 



The inclusion of the decuplet baryon fields in HB^PT requires additional operators in- 
volving powers of A/A x , as A is a chiral singlet. Thus all terms in the Lagrangian can 
be multiplied by arbitrary functions of A/A x without altering their symmetry properties 
under chiral transformations. To the order we are working, we can systematically include 
the contributions to the octet masses from these additional operators, by treating all con- 
stants in the calculation as arbitrary polynomial functions of A/A x , and expand them to 
the appropriate order. For example (Xm,I3m an d o"m must be expanded out to C(A 2 /A X ) 
while D, F and C must be expanded to 0(A/A x ), 



D ^ D 



A 

K 

A 



D 



1 + a 1 ; 

1 + 



A 

A v 



A 

A3 



+ o 



A 3 

x 



A 

A v 



A2\ 



0\ I 



(26) 



Additionally, the LECs are also implicit functions of /i, to cancel the scale dependence of 
observables arising from the loop integrals. 

In calculating the masses, the leading dependence upon m q arises from the terms in 
our Lagrangian, Eq. (fTHj) . with coefficients qzm,Pm, and o"m, which can be determined 
from experiment or lattice simulations. The 0(rriq^ 2 ) contributions arise from the one- loop 
diagrams shown in Fig. ^ which are formed from the operators in the Lagrangian with 



S 



(2) 

FIG. 2: One- loop graphs which give contributions to M B . Single, double and dashed lines 
correspond to octet baryons, decuplet baryons and mesons, respectively. The filled squares and 
filled circles denote axial couplings and insertions of given in Eq. (|15|) . The empty squares 
denote insertions of fixed 1/Mb operators from Eq. (J22J), and the empty circles denote two-baryon- 
axial couplings defined in Eq. H23|) . 

coefficients a, (3 and C. The 0(m 2 ) contributions arise from the one-loop diagrams shown 
in Fig. |2l from the tree-level contributions of the operators with coefficients, bf 1 , and from 
the NLO wave-function corrections. 

We find that the contributions to the nucleon mass are 5 



M 



(i) 



M 



N 

(3/2) 
N 



2m(a M + Pm + 2ajif) + 2m s a M , 

hD + Ffml + \{D-?>Ffml + ]- (5D 2 - 6DF + 9F 2 ) 

i >7TJ Z \_2 o '3 

C 2 (A 1 

H ^-^(^tt, A, n) + -T{m K , A, /i) 

n \ 6 6 



(27) 



K 



(28) 



consistent with previous calculations j23[ 36[. The contribution of this work is the NNLO 



5 We express a and (3 in terms of D and F, following the convention set by 1 



mass contribution to the octet baryons in 577(3), 



M 



(2) 



N 



-£(mmfi) 



(4vr/)3 
4 



l h A + l h A _ b A + 3b A + l h vA + 3 A _ 1 h A 3 A 

2 1 2 2 3 4 8 1 8 2 4 3 4 4 
27 



M b (4tt/) 2 



— (D + F) 2 + -C 2 
16 V ; 2 



-C{m v ,jj) 



(4tt/)» L6 1 6° 2 + 6° 3 + 64 + 24 61 + 24 ° 2 



mm 



24 



-C(m K ,ii) 



M B (A7rf) 2 16 
1 



(/J - 3F) ; 



(4tt/)» 
4 



b A + b A + ±b A + ^b vA + ^b vA + b vA 

1 r 3 (5jD 2_ 6jDF + g F 2 )+ 



M b (4tt/) 2 



m6(« M + /3 M + 2cr M ) 1 
+£{m 7n ii) h L[m^n) (A _ r 2 



(4tt/)s 



(4tt/) s 



-m(aM + /3m + 2cr M ) + -m« 



+£(m K ,n) 



(4tt/)= 



2(m + m s )(a M + /? M + 4a M ) + (5/J 2 - 6/JF + 9F 2 )M 
-^(/J-F^mW-^/J + SF) 2 ^} 



(i) 

AT 



+.7(771,,, A, /i) 



4C 2 



(4tt/)s 



+777, 



1 



+777 



+777 



* I (4tt/) 3 

4C 2 
^(4T/F 
1 



+ J(m K ,A,/i) 
+ 



C 2 



16 1 16 2 8 3 8 4 



M# } + Mi 1} 



(4tt/) 2 
1 



M$ + 



+777 



" I (4tt/) 3 
1 



Lm A + -ifef + J-6f + -b vA 
48 1 48 2 48 3 8 4 



K 



\h vA 4- -fr^ 4- -//^ 



+ 



(4^7f{^ (5jD2 ~~ 6DF + 9F2)M ^ ~ 3(jD ~ F)2M * 



M^l^ D2 - QDF ^ 



1 



-(/J + 3F) 2 M| 1} +C 2 
3 



M# } + M§ 



+ 



2 

777 



+ 



(4tt/) L 
mm 



(4tt/) 



6f + & 2 + &T + 26f + 26f + 26f + 4fef 

2 



6f +46? 



+ 



nr 



(4tt/) 



6f + 6f 



(29) 



The LO octet and decuplet masses given in the above expression can be found in AppendixlBl 
We have made use of several abbreviations for the functions entering from loop contributions, 
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namely: 



£(777,4,, //) 



777w log 



m 



F(m4,A,fi) = (ml -A 2 ) 



A 2 — ml lot 



A - J A 2 -ml + ie 



A + J A 2 -ml + ie 



( m \ 
A.og [J 



\ ( Tfl 

-gAmJlog 



2 / ' 



77 



2 A, /A 2 -m! lot 



A - 4 / A 2 - m 2 + ie 



A + ,/A 2 -m 2 +ie 



(30) 



It should be stressed that the above expansion, Eq. (|2*7|). and (j2SJ), is a quark mass 
expansion. The meson masses are replacements for the quark masses, given by Eq. (J5|). If an 
expansion in terms of the physical meson masses is desired, then to be consistent to £?(m 2 ) ~ 
0(mA), one needs to use the NLO expression for the meson masses in Eq. (J2T1) and the LO 
expression in Eq. (|2*5|) and Eq. (J29J1 . These results build on previous work [2y, Q, H3, EH . 
The results for the xPT calculation of the S, A and 5 baryons are located in Appendix IBl 



IV. PQxPT 

In PQQCD Lagrangian is 

9 

£= ^Q'W-motfQk. (31) 

j,k=i 

This differs from the SU (3) Lagrangian of QCD by the inclusion of six extra quarks; three 
bosonic ghost quarks, (u, d, s), and three fermionic sea quarks, (J,l,r). The nine quark 
fields transform in the fundamental representation of the graded SU{6\3) group. They are 
combined in the nine-component vector 




Here we have separated the quark field vector into the valence, sea and ghost sectors: 
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The quark fields obey the equal-time commutation relation 

Qf (x)Qf (y) - (-l)^Qf (y)Q?(x) = <T%5 3 (x - y), (34) 

where a,/3 are spin and i,j are flavor indices. Analogous graded equal-time commutation 
relations can be written for two Q's and two Q^s. The grading factors, 

/ 1 for fc = 1,2,3,4,5,6 . . 

% ~\0 for A; = 7,8,9 ' ^ 

take into account the different fermionic and bosonic statistics for the quark fields. In the 
isospin limit the quark mass matrix of £77(6 1 3) is given by 

rriQ = di&g(m,m,m s ,mj,mj,m r ,m,m,m s ). (36) 

Setting the ghost quark masses equal to the valence quark masses results in an exact can- 
cellation in the path integral between the valence quark determinant and the ghost quark 
determinant. The sea quark determinant is unaffected. Thus one has a way to vary the 
valence and sea quark masses independently. QCD is recovered in the limit rrij — * m and 
m T — > m s . 



Pseudo-Goldstone Mesons 



For massless quarks, the theory corresponding to the Lagrangian in Eq. (p?T]) has a graded 
SU (6\3) SU(6\3)r®U(1)v symmetry which is assumed to be spontaneously broken down 
to SU(6\3)v <S> U(l)v in analogy with QCD. The effective low-energy theory obtained by 
perturbing about the physical vacuum state of PQQCD is PQxPT. The result is 80 pseudo- 
Goldstone mesons whose dynamics can be described at LO in the chiral expansion by the 
Lagrangian 



£ = t- s tr (O^d^E) + Astr (m Q E + m^E*) + a*0"$ o fy$o 



where 



exp 



2z$ 
T 



M X t 
X M 



(37) 



(38) 



The operation str( ) in Eq. (J3*7|) is the supertrace over flavor indices. The quantities a$ and 
mo are non- vanishing in the chiral limit. M and M are matrices containing bosonic mesons 
while x an d X are matrices containing fermionic mesons: 



M 



X 




J° 

J- 

J- 



L+ R+\ 
L° R° 

L° R° a 



Y, 



Y, 



jr 



J l 

m 

Y° 

1 lr 



■ jr 
Vr ) 





( Vu 






M = 




Vd 


k° 






~ 








K 


Vs J 



Xk+ Xj° Xl+ Xr+ 
Xko Xj- Xlo xrp 
Xvs Xj- Xlo xro 



(39) 
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The upper 3x3 block of M is the usual octet of pseudo-scalar mesons and the remaining 
components are mesons formed with one or two sea quarks. 

The flavor singlet field is defined to be $o — str($)/v / 6- PQQCD has a strong axial 
anomaly, U(1)a, and therefore the mass of the singlet field, m can be taken to be the 
order of the chiral symmetry breaking scale, mo — * A x [2l|. In this limit, the rj two-point 
correlation functions deviate from the simple, single pole form. The r\ a r\ b propagator with 
2 + 1 sea-quarks and a,b = u, d, s, j, I, r, u, d, s is found to be at leading order: 



ie a Sab i e a e b (q 2 - m%) (q 2 - m 2 rr ) 



VaVb y 2 _ m 2^ _j_ i e 3 ^2 _ m 2^ _|_ ^ ^2 _ m 2^ _|_ ^ (g2 _ m 2^ _)_ 

where 



(40) 



e a = (41) 

The mass, m xy , is the mass of a meson composed of (anti)-quarks of flavor x and y, while 
the mass, mx is defined as m 2 x = ~ (m|- + 2m^ r ). This can be compactly written as 



where 



'VaVb 



eJabPa + e a e b H ab (P , P b , P x ) , (42) 



p _ p p 

a 9 9 i • ' O 9 9 i • ' X 9 2 i • 

g z — mi + ze q z — + le q z — + le 



n ab (A,B,C) 



1 

3 



{m% - m 2 J (m 2 . r - mjj ^ _ (m^. - m 2 J (m 2 „ - m 2 J 



ru- 



in: 



)( 



Hi 



'/a 



m 



x 



) ( 



'/a 



m 



m 



)( 



m 



x 



-) 



-m 2 a ) (m 2 x -m 2 b ) 



C 



[m 



x 



m 



in, 



)( 



m 



'ib 



m 



x 



B 



(43) 



B. Baryons 

In PQxPT the baryons are composed of three quarks, QiQjQk, where % — k can be 
valence, sea or ghost quarks. One decomposes the irreducible representations of SU (6\3)v 
into irreducible representations of SU (3) m / ® SU(3) sea ® SU(3) g h ost ®U(1). The method for 
including the octet and decuplet baryons into PQxPT is to use the interpolating field [l3l.l32^: 

B] jk ~ (QT a Qf b Ql' c - QT a Q]' c Qk b ) ^(^W (44) 

We require that B^ k = Bij k , defined in Eq. (fTTj). when the indices, i, j, k are restricted to 1—3. 
Thus the octet baryons are contained as an (8, 1, 1) of SU (3) va i<g)SU (3) sea ®SU (3) g h os t®U(l) 
in the 240 representation. In addition to the conventional octet baryons composed of valence 
quarks, Bij k also contains baryon fields composed of sea and ghost quarks. In this paper we 
only need the baryons which contain at most one sea or ghost quark, and these states have 
been explicitly constructed in Under the interchange of flavor indices, one finds [32^ : 

B ijk = (-) 1+ ^B lk , and B ijk + (-) 1+ ^B jik + (-y+ViVi+ViVk+ViV*^.. = . (45) 
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Similarly, one can construct the spin-| decuplet baryons which are embedded in the 138, 
and have an interpolating field 

%■£ ~ (Qt' a Qf b Ql' c + Qi' b Q]' c Qk a + QT c QT a Q" k b ) e abc (cy ^ • (46) 

We require that %j k = Tij k , when the indices i, j, k are restricted to 1 — 3. Under SU(3) va i (g> 
SU (3) sea <E> SU(3) g host ®U(1) they transform as a (10, 1, 1). In addition to the conventional 
decuplet resonances composed of valence quarks, %j k contains fields with sea and ghost 
quarks, which have also been constructed in [lj|. 

Under the interchange of flavor indices, one finds that 

%jk = {-) 1+mVj T jik = (-) 1+ ^ k T ikj . (47) 

Under SU(6\3)v, both Bij k and %j k transform as 32] 

B ijh -> (-)^ + ^+^ + ^^+^U/U j m U k n Bi mn . (48) 



C. Baryon-Meson Lagrange Density 

The utility of using 3-index baryon tensors in the QCD calculation is fully realized when 
one extends the %PT Lagrangian to PQxPT. It has recently been shown that in the meson 
sector there are extra operators in PQxPT, that do not arise in %PT ^3]. The Cayley- 
Hamilton identity for traceless 3x3 matrices, allows one to reduce the number of operators 
in the meson sector of xPT. There is no such identity for the supertraceless matrices of 
S77(6|3), allowing for extra operators. Use of the SU(3) matrix of baryon fields, Eq. (JTUJ), 
raises similar concerns at higher orders in the expansion of the Lagrangian, when extending to 
PQxPT. However, when using three-index baryon tensors, the flavor contractions invol ving 
the baryon fields can not be expressed as a trace, as with the more familiar form ji^, I24T 
I25L I26L I341 . 136^ . Therefore, there is no danger of inadvertently neglecting operators when 
extending from £77(3) HB^PT to £[/(6|3) PQHBxPT. Furthermore, when employing the 
three-index notation for the baryon fields the coefficients of the HBvPT operators are the 
same as the coefficients of the corresponding PQHBxPT operators j39|. This is useful as for 
the foreseeable future, reliable calculations of baryonic matrix elements will be performed 
with PQQCD lattice simulations. 

To write down the PQ%PT Lagrangian, we must also include the appropriate grading 
factors. The flavor contractions are now defined as lli 



{By b) — W Ml vi B Mk , (T r % ) = T^ kji r mjh 

(byyb)= T' kii y { 1 B Mk , (fry %) = r rti y ;i 1 r mik 

(BY BY) = (-Y^+VjKvk+m)^^ r £ y k l B Pjijl , 

(B Y Y» T M ) = T' m Yi (Y»\ l r mjk . (49) 
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The leading order PQ Lagrangian is given by 

££g = (BivDB) + 2a M (BBM+) + 2(3 M (BM+B) + 2o M (BB) str(M+) 
-(7"[ivD- A]TJ + 2 7M (T^ + r M ) - 2a M (r M Tj str(X + ) 
+2a (JBSPBAp) + 2/3 (BS^B) + 2H(T v S f *A t X) 

+ ^-C[(TA u B) + (BA U T U )] . (50) 
The fixed coefficient and higher dimensional operators Lagrangians are given by 

£*? = - (g-5l B ) - a (b^Bv .A\ +a (b^Ibv ■ A 
we \ 2Mb J y Mb J \ M B 

and 6 

£%? = ^yyj b A B kji {A-A\ n B njk + 6^ . A) k n B ijn 

+ b A(_y l (r,i+v n )B kj \AJ i l {A% n B lnk + b A B kji B ijk Tr (A ■ A) 

+ b vA B kjl (v-Av A); B njk + bf (-)(.vi+%)(v k +v n )B kji (y • Av • A) k n B ijn 

+ bl A (-)^+^B kj \vA) i l (vA) j n B lnk + bl A B kji B ijk Tr:(vAvA) 

+ b A 1 I B k3 \M + M + ) t n B n , k + b^(-)^ + ^ + ^m kj \M+M + ) k n B ijn 

+ b f (-)m(n 3+ v n )B kj \ M+ ) t l (M + ) J n B lnk + fof B kji B ijk str (M + M + ) 

+ bf B kj \M + ) i n B njk str (M+) + 6f B^By* str (A4+) str (M+) 

+ be {-) { ^ ){Mn) B k3 \M+) k n B ijn str (M+) } . (52) 

All the coefficients appearing in Eq. (|5UJ1 . Eq. (|51|) and Eq. (|52j) . have the same numerical 
values as in x?T. The situation is similar to the case considered in Section III Bl 

All the operators in the Lagrangian can be multiplied by arbitrary functions of A/A x . We 
include these effects by treating all the coefficients in the Lagrangian as arbitrary polynomial 
functions of A/A x , and expand out to the appropriate order (see Section All the 

coefficients appearing in this expansion also have the same numerical values as in the %PT 
case. 



V. BARYON MASSES IN PQxPT 

The chiral expansion of the octet baryon masses in PQxPT has the same form as in ;\PT. 

M m = M 0*) - M« (/x) - M| /2) (/,) - Mg> (/.) + ... (53) 

6 It is interesting to note operators of the form (B[A^, A v \ S v \B) that vanish in xPT don't vanish in 
PQxPT due to grading factors. However, they don't contribute to the octet baryon masses. 
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FIG. 3: In addition to the one- loop diagrams in Fig. ED Mg' ' also receives contributions from the 
singlet (hairpins) in PQxPT. Single and double lines correspond to 240-baryons and 138-baryons 
respectively. The crossed dashed line denotes a hairpin propagator. The filled squares denote the 
axial coupling given in Eq. (|5Ujl. 






FIG. 4: In addition to the one-loop diagrams in 

Fig.0 M|. also receives contributions from the 
singlet (hairpins) in PQ^PT. Single and double lines correspond to 240-baryons and 138-baryons 
respectively. The crossed dashed line denotes a hairpin insertion. Filled squares denote the axial 
coupling given in Eq. (|5Uj). Filled circles denote a coupling to Ai+ given in Eq. ()5U[). Empty squares 
and empty circles denote insertions of fixed 1/Mb operators given in Eq. 1)51(1 . and insertions of a 
two-baryon-two-axial coupling given in Eq. 1)52)1 respectively. 



However, we must also include hairpin contributions from the eta field propagators. See 
Fig. El and Fig.H 

We find that the mass expansion of the nucleon is: 



M 



(i) 



M 



N 

(3/2) 
N 



2m(a M + Pm) + 2a M {2rrij + m r ) 

— AD{F - l -D)ml + ^(hD 2 - 6DF + 9F 2 )(2m% + m 3 n 

2C 2 

+(D - 3F) 2 M 3 (m nj m,) + — T{m^ A, /i) 

2C 2 C 2 
+ ^—F{m ju , A, fi) + —T{m ru , A, /i) 

07T 6 



(54) 



(55) 



consistent with [13]. The contribution of this work is the NNLO mass contribution to the 
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octet baryon masses. 
-A4 2) = C{m^Li) 



(4tt/)3 



°3 + g °3 



M b (4tt/) 2 



3D(3F - £>) 5C 2 
2 + ~T 



(4tt/) j 



fc A + b A + b A + \b\ A + \bl A + \bf 



+ 



9(D - 3F) : 
8M b (4tt/) 2 



-£(m ju ,Li) 



1 



(4tt/)3 
4 



-£(m ru ,Li) 



M B (4nf) 2 
1 



i 2 2 1 2 2 
1 f3(5D 2 - QDF + 9F 2 ) 5C 2 



(4tt/)s 
4 



(47T/)* 



M b (4tt/) s 
1 . 
4 



b A + b A + \b\ A + \bl A 



-(5D 2 -6DF-9F 2 ) + -C 2 

8 8 



h A 4- -h vA 
o 4 -I- o 4 



2(mj + m r )C(mj r , fj) + m 
1 



4<r 



A4 



(4vr/) 2 
+C(m ju ,ii) 



(4tt/) s 



+C(m ru ,/j,) 



(4tt/) s 



2C 2 



4£(m % , //) + 4£(m jr , /i) + C(m Vr , fx) 
+2C(m Vj , m^. , /i) + £(m^ r , , //) 

2£(m % , /i) + £(m % , m % . , //) 

4(m + m 3 )(a M 4- /3m) 4- 2(5D 2 - 6DF + 9F 2 )M^ 1} 

-9(D-F) 2 M«-(L>4-3F) 2 M«} 
2(m + m r )(a M 4- /3 M ) 4- (5D 2 - 6DF + 9F 2 )M { ^ 



-\{D-F) 2 M^- l -{D + m 2 M^ 



+J(m ru ,Li) 



(4tt/) 2 

_C 2 _ 

(47T/) 5 



+ Mi 1} 



2C 2 



M^ 1} + Mg 



+ 



2(4tt/) j 



mi. + mi r + -m* 

<h J' A 'It 



1 b VA 



+ 



7r U47r/) 3 16 M b (4tt/) 



5D(3F-D) 9C 2 



4-m r 



+ 5(5L> 2 - 6L>F + 9F 2 ) + 9C 2 

4(4tt/) 3 8M b (4tt/) 2 
6f + 6^ 5(5D 2 - 6DF + 9F 2 ) + 9C 2 



8(4tt/)= 



16M B (47r/) ; 



+ cont. 
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b vA + b vA + bl A 15(D - 3F)' 



lQM B (A7rfy 



h vA 



K4tt/)3 
2 



+m 



+m 



2 



2.M {m v .,m m ) + M (m Vr ,m Vr ) 



(4tt/) s 



;(5D 2 - 6DF + 9F 2 )M^ - 3(D - F) 2 m£ ] 



1 



-p + SF^M^+C 2 



(1) 



■(4tt/)= 



-(5L> 2 - 6DF + 9F' 2 )M^ ] - 3(D - F) 2 M^; 
3 



(i) 



i 



-(D + SF^M^+C 2 



(i) 



+ 



2 

m 



M 

i 



6f + 6f) + 



(4tt/) 



mmj -(26f + 26f ) + 7^(26f + 46f ) 



(4tt/) 



(47T/) 



(6f + b 



where the functions not already defined in Eq. (|30|) are 



f), 



Hw (C(m<f„ fi), L{my, /x), £(m x , /x) j 
7^ ( £(m^, /x), £(rry, /x), £(m X) A*) 



(56) 



(57) 



We stress that these mass expressions, Eq. (|54 J1 - (|56J1 . are quark mass expansions. If a meson 
mass expansion is desired, one can express these quark masses in terms of the lattice meson 
masses. To do this consistently, one must replace the quark masses in Eq. (I54j) with the 
NLO relation to the PQxPT meson masses, while the quark masses in Eq. (|B3j) and Eq. (JSTIj) 
can be replaced by the LO relation to the PQ meson masses. It is easy to show that these 
expressions reduce to the expressions, Eq. fl2Tj) - (}2T?j) . in the limit that rrij — > m and 
m r — > m s . The results for the NNLO mass calculations for the remaining octet baryons can 
be found in Appendix O 



VI. CONCLUSIONS 

We have calculated the (9(m 2 ) contribution to the masses of the octet baryons in the 
isospin limit of and PQxPT, keeping the decuplet baryons as dynamical intermediate 
states. Working to this order in the chiral expansion introduces a large number of LEC's. In 
order for the calculation to have any predictive power, these LEC's must be fit from various 
experiments and or lattice results. Lattice calculations will eventually allow first principles 
determination of these constants, and thus predictions of QCD observables. However, for the 
time being, to make rigorous predictions of the baryon mass spectra arising from QCD, one 
needs to perform PQQCD lattice simulations and extrapolate these results using PQ^PT. 
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To date, there are limited PQQCD lattice calculations of the octet baryon masses. We 
hope this work will help to further interest an execution of these calculations. This is the 
first calculation to 0(m 2 ) in the baryon sector using PQ^PT. Calculations of this order 
are necessary for reducing the uncertainty in the chiral extrapolations and systematically 
studying the chiral expansion of both PQxPT and x?T. 

Note Added. While this work was being completed, an independent analysis by Frink and 



Meissner appeared 41]. Their work considers the octet baryon masses in x?T only. 
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APPENDIX A: CHIRAL LOGS AND OTHER FUNCTIONS 

For the readers convenience, here we list all the functions arising in the calculation of the 
octet baryon masses, including the previously defined ones (which have been left unnum- 
bered). 



£(m0,/x) = ml log 



Cim^n) = m log — 



mz 



mz 



T{m^ A, n) = (mj - A 2 ) 



A 2 — m 2 log 



A 



i?) 



— A///,.loR I — j 



A 



A 2 - ml + ie \ / m 2 
-Alog^ 



A + y'A 2 - m 2 + ie 

2 



+2A./A 2 - m 2 , log 



A 2 — m\ + ie 



A + JA 2 -m 2 +ie 
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M n {m^m^) = 


Ti.^ 1 


(m^,m^,m x ) 


£(777,4,, 77V ) = 




( £(mfa ft), £(171^,(1), C(m x , fi)j 


Lim^vn^,) = 




[C(mfa fi), C(m^, (jl), C(mx,/j)\ 


T{jn$,mtf) = 




A, fi), T[m^, A, fi),T{m X i A, fi)j 


Jim^m^,) = 


T~C<f><f>i 


\J{m^ A, /jl), J(m<y, A, //), J(m x , A, //) j 


And H a b(A, B,C) is denned in Eq. 


(USD. 



(Al) 
(A2) 



APPENDIX B: OCTET BARYON MASSES IN %PT 

The mass of the i th baryon in the chiral expansion is 



M Bi = M (/i) - M« {fi) - M% A > (fi) - M% (/x) + . . . (Bl) 

The LO and NLO mass corrections to the octet baryons are listed for completeness. The 
LO mass corrections are listed in Tabled The NLO mass corrections are given by, 



r(3/2) 



r(2) 



M 



(3/2) 



(4tt/)' 



E 



(G BB ^m\ + C 2 {C TB ^fT{ m!j) , A, y) 



The sum on <p is over loop mesons with mass m^. The C BB $ and are Clebsch- Gordon 

coefficients between the Axial field and two octet baryons or an octet and decuplet baryon. 
The sum of these coefficients are listed in Table |TI] 

Results of this work are the NNLO mass corrections in to the octet baryons given 
by the expression, 

_l_>^r„„ 1 
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Si 
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1 3 



n A _i rr 



3C(m^,fi) + 



(47r/) 2 8MB E{( c '^) 

+C ,2 (Cr J B <} i) 2 
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Cim^fi) + -mj 
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C{m^,fi) + - 



(4tt/)= 



J(m^,A,fi) + 



J (m^, A, /i) +m 2 



(4tt/) 



.1714171 j 
3 1 J 



(B2) 
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The coefficients for the masses are listed in Table ITT1 through Table rvTTl 

TABLE I: The LO Octet and Decuplet Baryon Masses in xPT. 





Octet Baryons 

<> 




Decuplet Baryons 
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2m(a M + Pm + 2g m ) + 2m s a M 


Mi" 


2m(7M - 2o ; j W -) - 2m s a M 




m (|«m + IPm + 4<T M ) 
+m s {\a M + |/?m + 2ctm) 


M« 


|(2m + m s )(7 M - 3ct m ) 
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m(a M + 2f3 M + 4ctm) + rn s (a M + 2a M ) 


M« 


|m(7M - 6a M ) + §m s (27 M - 3a A f) 




m (\a M + |/3m + 4fJA/) 
+m s (|q m + |/5m + 2<tm) 







TABLE II: The octet-octet-axial and decuplet-octet-axial coupling coefficients in xPT. The coef- 
ficients have been grouped according to loop mesons with mass for each octet baryon. 
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TABLE III: The coefficients of 2-octet 2-axial couplings in xPT. The coefficients have been grouped 
according to loop mesons with mass for each octet baryon. 





7T 


Yl(p ^BB4><i> an< ^ S</> ^ V BB<i><\> 
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TABLE IV: The C BB ^ coefficients in xPT. The coefficients have been grouped according to loop 
mesons with mass for each octet baryon. 





riM 
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TT 7] K 


TV 
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TABLE V: The coefficients for the internal octet mass insertion and the wavefunction correction 
in xPT. The coefficients have been grouped according to loop mesons with mass for each octet 
baryon. 
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TABLE VI: The coefficients for an internal decuplet mass insertion and the wavefunction correction 
in xPT. The coefficients have been grouped according to loop mesons with mass for each octet 
baryon. 
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TABLE VII: The Cb 4 coefficients in xPT. The coefficients are grouped by products of quark 



masses for each octet baryon. 
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APPENDIX C: OCTET BARYON MASSES IN PQxPT 



The mass of the i th baryon in the chiral expansion of PQHBxPT has the same form as 
in HBxPT. 



M Bi = Mo 0*) - M« 0*) - M%'> (/*) ^ M^ 2 / (/i) + 



(3/2) 



r(2) 



(CI) 



The LO and NLO mass corrections to the octet baryons in PQxPT are listed here for 
completeness, and are consistent with The LO mass corrections for the 240- and 

138-baryons are listed in Table IViLTl Table IXl and Table HXl The NLO mass corrections are, 
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(3/2) 



(47T/) 2 
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(4tt/) 



2E 



(Cb^) 2 ™J + G 2 (^) 2 ^K, A, ,x) 



(C2) 



The sum on is over loop mesons with mass m^. The and are PQHBxPT 

Clebsch- Gordon coefficients between the axial field and the baryon fields. The sums on <p, <// 
are over contributions from the hairpin eta field interactions. The sum is over pairs of eta 
fields, with no double counting. 
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The main results of this work are the NNLO PQHB%PT octet mass calculations. 
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The coefficients for the NLO and NNLO octet baryon masses are listed in Table IXTl through 
Table IXXII It is straight forward to show that in the limit rrij — > m and m r — > m s , these 
expressions reduce to the xPT expressions listed in Appendix |Bj 
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TABLE VIII: The LO 240 (B) and 138 (T) baryon masses composed only of valence quarks in 
PQ X PT. 





Octet Baryons 
, ,(11 

M B 




Decuplet Baryons 
, ,(il 




2m{aM + Pm) + ^vrijOM + 2m r <j M 


Mi" 


2m~fM — ^rrijaM — 2m r aM 


*?' 


m (\a M + IPm) + m s (\a M + |/?m) 
+4mj(JM + 1m r OM 


M<V 


ffnjM + \m s lM ~ ^mjWM - 2m r a M 


«1" 


m(aM + Wm) + m s aM + ^m^ou + 2m r a M 


m<:» 


ImjM + |"t-s7m - kmjOM - 2m r a M 


a4" 


m {\a M + f /3m) + m s {\a M + |/?m) 
+Amj(JM + 2m r <JM 


<» 


2to s 7m — 4m J o r M — 2m r aM 



TABLE IX: The LO masses of the 138 (T) baryons composed of 2 valence quarks and 1 sea or 
ghost quark in PQxPT. 
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TABLE X: The LO masses of the 240 (B) baryons composed 2 valence quarks and 1 sea or ghost 
quark in PQxPT. 
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— 1. 1 / Udl V Ullo 

with one ghost quark 




m (|a M + §/3m) + 2m r a M 
+mj [\a M + f /?m + 4a M ) 


Z-'u 




<> 


to (|a M + |/3 M ) + AmjaM 
+m r (\a M + §/3m + 2cr M ) 






<> 


m(a.M + 213m) + mj(aM + ^<?m) + 2m r aM 


/l \ 

M A 






m(a M + 2(3m) + ^rrijUM + m r (aM + 2<jm) 




Mi" 




(to + to s ) (|a M + §/?m) + 2m r a M 
+mj (\a M + \(5m + 4<7 M ) 


\ 


to (|«m + §/3m) + 4mj(7M 
+m s (§a M + |/3m) + 2m r a M 


n ) 


(to + to s ) (|a M + §/3m) + 4toj0m 
+m r (\a M + §/3m + 2cr M ) 


JW g(6) 
"a 


to (|a M + |/3m) + 4toj0m 
+m s (la M + §/3m) + 2m r a M 


"% 


±(to + m s )(a M + Wm) + mj(a M + 4<tm) 
+2m r <TM 


'u 


fn (|a M + /?m) + kmjOM 
+m s (\a M + /3m) + 2m r a M 


M (1) 

= (3) 


i(m + m s )(a M + 2/3m) + 4toj<jm 
+m r (aM + 2a m) 


M {1) 

= (3) 


fn {\a M + /3m) + 4:mja M 
+m s (|«m + /3m) + 2m r a M 




m s (|a M + §/3m) + 2m r a M 
+mj {\a M + \$m + 4ctm) 


< 




<> 


to s (|«M + §/3m) + kmjVM 
+m r {\oim + §/3m + 2gm) 


<> 


Amj(jM + 2m s (a.M + Pm) + 2m r UM 



TABLE XI: The 240-240-Axial coupling coefficients in PQxPT. The coefficients are grouped by 
the loop mesons with mass for each of the octet baryons. 
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TABLE XII: The coefficients for meson loops containing hairpin insertions with an internal 240- 
baryon. 
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TABLE XIII: The 240- 138- Axial coupling coefficients and the coefficients for meson loops con- 
taining hairpin insertions with an internal 138-baryon in PQ%PT. The coefficients are grouped 
by the loop mesons with mass and by the participating 77 fields in the hairpin interaction, 
respectively. 
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TABLE XIV: The and C^'J? coefficients in PQxPT. The coefficients are grouped by the 

loop mesons with mass m^. 
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TABLE XV: The C^Z, and C B aJ> coefficients in PQxPT. The coefficients are grouped by the 



77 fields participating in the hairpin interaction. 
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TABLE XVI: The C B ^^ coefficients in PQxPT. The coefficients are grouped by the loop meson 
with mass m^. 
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TABLE XVII: The C ' BB<jxj>' coefficients in PQxPT. The coefficients are grouped by the rj fields 
participating in the hairpin interaction. 
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TABLE XVIII: The coefficients for internal 138-baryon mass insertions and the wavefunction 
corrections. The coefficients are grouped by the loop mesons with mass for each octet baryon. 
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TABLE XIX: The coefficients for hairpin interactions with internal 138-baryon mass insertions 
and wavefunction corrections. The coefficients are grouped according to the r] fields participating 
in the hairpin interaction. 
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PO 

TABLE XX: The C B coefficients in PQxPT. The coefficients are grouped by products of quark 
masses. 
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TABLE XXI: The coefficients for internal 240-baryon mass insertion and the wavefunction corrections. The coefficients are grouped 



according to the loop mesons with mass for each octet baryon. 
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